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Abstract 

We study the stability constraints on the parameter space of a triplet extension of MSSM. 
Existence of unbounded from below directions in the potential can spoil successful Electroweak 
(EW) symmetry breaking by making the corresponding minimum unstable, and hence the model 
should be free from those directions. Avoiding those directions restricts the parameter space of 
the model. We derive four stability constraints, of which only three independent from each other. 
After scanning the model’s parameter space for phenomenologically viable data points, we impose 
the stability constraints and find that only about a quarter of the data points features a stable 
EW minimum. At those data points featuring stability, y and the up Higgs soft mass turn out to 
be smaller than about a TeV in absolute value, which make the mass of the lightest chargino and 
neutralino smaller than about 700 GeV. Two relevant phenomenological consequences of lifting the 
unbounded from below directions are that the lightest Higgs boson decay rate to diphoton predicted 
by the triplet extension of MSSM generally features larger deviations from MSSM and fine tuning 
is actually higher, that what each of the two would be without imposing stability constraints. 
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I. INTRODUCTION 


The Triplet Extended Supersymmetric Standard Model (TESSM) was introduced mainly 
to enhance the tree level Higgs boson mass while satisfying the top-quark mass bound [1, 2]. 
The authors have shown there that Electroweak (EW) symmetry breaking is successfully 
realized in this model, and studied constraints on the parameter space which must be satisfied 
for the potential to be stable. The lightest Higgs boson mass has been calculated in [1-3] 
for the Minimal Supersymmetric Standard Model (MSSM) with extra scalar triplet chiral 
superfields with hypercharge, respectively, Y = 0, ±1. The one loop correction to the Higgs 
boson mass was calculated for MSSM with a Y = 0 scalar triplet also in [1], The authors 
there pointed out that a light Higgs boson mass of 0(100) GeV can be generated already 
at tree level, if the triplet coupling to a pair of Higgs bosons is large and comparable to the 
top-quark Yukawa coupling. More recently it has been shown that the triplet charged states 
in TESSM can comfortably enhance the diphoton decay rate of the Higgs boson to match 
the value observed at LHC [5, 6]. Further phenomenological studies of TESSM explored 
neutrino mass generation and leptogenesis [7], charged Higgs production at colliders [8], 
spontaneous CP violation [9], etc. 

In this paper, we are interested in studying the stability of the EW minimum of the 
TESSM scalar potential. If the EW minimum is not a global minimum, correct EW symme¬ 
try breaking is not realized and its viability spoiled. 1 It is therefore important to determine 
the constraints on the parameter space that ensure that the EW minimum is stable. For 
MSSM, there are a few directions possible along which the potential becomes Unbounded 
From Below (UFB), as it has been shown in [10]. The authors there also discussed the 
Charge and Colour Breaking (CCB) minima for MSSM. Other relevant studies on unstable 
and metastable minima of the supersymmetric scalar potential can be found in [ 1-16]. The 
problems associated with this kind of minima were addressed in [17]. The study of unreal¬ 
istic vacua or CCB minima in different supersymmetric models have already been discussed 
in [18] for MSSM with neutrino mass operators, in [19, 20] for NMSSM, in [21] for i/SSM 
etc. 

For TESSM, conditions for the stability of the potential have been derived in [ ]. Indeed, 
as we have pointed out above, there may be a few unaccounted for UFB directions in the 

1 We assume throughout this paper the age of the Universe to be infinite. 
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field space. Hence we would like to perform a full analysis of the UFB directions of the 


TESSM potential. We shall show that stability constraints that lift these directions allow 


one to constrain severely the parameter space of TESSM, with observable consequences for 
the mass spectrum and TESSM phenomenology. 

This paper is structured as follows: In Section II we introduce the model, define the EW 
minimum of the potential, and discuss the present theoretical constraints. Next, we present 
the main results of this paper, which are the UFB directions of the TESSM potential and the 
corresponding stability constraints on the EW mininum, in Section III. In the subsequent 
Section IV we impose the stability constraints on a large set of data points, satisfying the 
present experimental constraints, and show how various relevant quantities are affected by 
the new constraints. Finally we draw the conclusions in Section V. 

II. THE MODEL 


The held content of TESSM is equal to that of MSSM extended by a Y — 0 SU(2) triplet 
chiral superfield, whose scalar component can be written in matrix form as 



( 1 ) 


The renormalizable superpontential of TESSM includes only two extra terms as compared 
to MSSM 2 , given that the cubic triplet term is zero: 


Wtessm — ^tTt(TT) + nc>Hd-H u + A Hd-TH U + y t UH u -Q — ybDHd-Q — y T EHd-L , (2) 


where represents a contraction with the Levi-Civita symbol e iv with e 12 = —1, and a 
hatted letter denotes a superfield. For later convenience we write explicitly the potential D 
and (quartic) F terms derived from the superpotential above: 



( 3 ) 


( 4 ) 


2 To simplify the phenomenology of the model we neglect the first and second generation Yukawas, as well 


as R-parity breaking terms. 
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where T are group generators with gauge indices a and b , and cjfj runs over the scalar 
components of all the TESSM chiral superfields. 3 Also, g 3 , g l, and r/y are the SU(3)c, 
SU(2) l , and the f/(l)y gauge couplings, respectively. The full potential is then given by: 

V = V D + V F + V s , (5) 


where Vs represents the soft supersymmetry breaking terms corresponding to the superpo¬ 
tential in Eq. (2), as well as the soft squared masses associated to each scalar field: 


Vs 


HtBtTt(TT) + gE>Bc,Hd-H u + A AtH,i-TH u + y t A t t* R H u -QL + h.c. 


+rc4Tr(T f T) + m 2 Hu \H 11 \ S + m\ d \H d \ 1 + ... . 


( 6 ) 


To break correctly the EW symmetry, SU(2) L xU(l)y —* U(1 )em, we assign nonzero vevs 
only to the neutral scalar components and impose the usual minimization conditions on the 
potential in Eq. (5) 

H u = ( a u + ib u ) , H% = -j=(a d + ib d ) , T° = -J= (a T + ib T ) ; 

^aj^lvev 0 , (^z) Vi ? ^ Zl^ d,T , (7) 


where a,- and A are both defined to be real. The conditions above allow one to determine 


three of the Lagrangian free parameters in terms of the other ones: 

2 2 9y “f / 2 2 \ I D Vd A / 2 I 2 \ , \ 

-5-K ~ v d)+ B d9d— - ~r (v d + Wj + A v T 


m H„ = -»D + 


2 , 9y + 9l / ,2 




1) + B D 9D ~ - ~r (vl + v%) + Xv T 

”d ^ 


A 2 


v.. — V. 


, At \ v d 


At \ 

“W 1“ P-T ) - 

2 / v d 


9d 


A 


riln 


- - — (y 2 d + u 2 ) - 2/r T (i? T + 2/i T ) + A 


Ad- 


2vt 


Ai 


- \ — + ^ t 


VdVu 

Vt 


( 8 ) 


By plugging Eqs. (8) in Eq. (5) one derives the expression for the potential at the EW 
minimum: 

" 2 1 - A 2 r o o o/o om AU- 


^EW = - 9Y ^ 2 9L i V d ~ V u) 2 ~Y + + A 2 J] 


[WW (A t + 2g T ) - (v 2 d + u 2 ) /r D ] 

(9) 


The EW potential given above receives relevant corrections at one loop, which can in prin¬ 
ciple be minimized by choosing a suitable renormalization scale close to the heaviest masses 


3 In Eqs. (2,3) as usual we dropped all the gauge and family indices that can be contracted, implied 
summation over repeated indices, and denoted with a tilde the scalar superpartner of a SM fermion field. 
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in the particle spectrum. Contrarily to MSSM, though, it is not possible to solve the EW 
vevs in terms of the couplings and dimensional parameters of TESSM, and therefore in 
principle one should run the vevs by using the Callan Symanzik equation for the effective 
potential of a softly broken supersymmetric theory [ ]. Analogously to the NMSSM case 

[18, 19], though, we choose to simplify our analysis by evaluating Vew at the EW scale 
v w = 246 GeV. In Section IV we find that the stability constraints obtained by comparing 
unrealistic vacua with Uew(w) are generally conservative. 

The first stability condition for successful EW symmetry breaking is obtained by requiring 
the trivial vacuum at the origin to be unstable. By taking all the vevs to be zero, the 
requirement that one of the eigenvalues of the neutral scalar squared mass matrix be negative 
is equivalent to imposing the condition 


B 2 d > III 


m 


H d 


i4j 


+ i 


m 


Hu 


9d 


+ 1 


( 10 ) 


m ? 0 < mi ( cos 2 2/3 H— 75 -——v sin 2 2/3^ , tan (3 
V 9y + 9l J 


( 11 ) 


When the condition above is satisfied, one can derive an important bound on the mass of 
the lightest neutral Higgs [1, 2]: 

\ 9 \ 

Vy 
Vd 

The result in Eq. (11) shows the main advantage and motivation of TESSM over MSSM: for 
tan/3 close to one and a large A coupling it is in principle possible in TESSM to generate 
the experimentally measured light Higgs mass already at tree-level [ 1], which would imply 
no or negligible Fine Tuning (FT) of the model. 

Even when the constraint in Eq. (11) is satisfied, for any given values of the free param¬ 
eters there can be colour and electromagnetic charge breaking minima that are deeper than 
the EW one or even unbounded from below (UFB) directions in the potential: given that 
the latter possibility gives the tightest constraints on the MSSM parameter space [10], by 
analogy in this work we focus our attention on UFB directions in the TESSM potential and 
the corresponding stability constraints, which we derive in the next section. 


III. UNBOUNDED FROM BELOW DIRECTIONS 

In softly broken supersymmetric models the potential is generally stable, given that the 
quartic terms are generated by the superpotential as well as by the gauge interactions (D- 
terrns) and the supersymmetric tree-level potential is semidehnite positive. If the quartic 


5 





terms cancel, though, the soft mass squared terms can eventually drive the potential to 
negative infinite values. Our aim in this section is to perform first a complete analysis of 
the possible UFB directions in the TESSM tree-level potential, Eq. (5), and then to derive 
the corresponding stability constraints on the parameter space. 

In general to find the deepest UFB direction of a supersymmetric theory in an N field 
subspace, one solves the minimization conditions with respect to N — 1 fields, and then 
substitutes the solutions in the potential which turns out not to have quartic terms. The 
combinations of vevs that can do the trick are those that can cancel separately each D 
and F quartic terms, Eqs (3,4), given that each of these terms is semidehnite positive. A 
straightforward way to cancel the SU(3) D-terms is to take the vevs of the left-handed (LH) 
and right-handed (RH) squarks to be the same. Analogously, we take the positive and 
negative charged triplet component vevs equal to each other, which cancels their SU(2) D- 
terms, and we do the same for the charged doublet components of H u and Hj, which cancels 
their U(l)y D-terms. To simplify our analysis we define also the LH and RH charged 
slepton vevs to be equal to each other. Finally, the cubic F terms, being supersymmetric, 
should be zero when the quartic terms involving the same fields are zero. We avoid cubic F 
terms from the outset by choosing sneutrino and slepton nonzero vevs to belong to different 
generations. In such a case, the surviving terms in the potential involving a stau vev can be 
readily obtained by simply relabeling those involving a sbottom vev, 4 and for this reason in 
the present analysis we simply neglect the stau vev and derive results involving it directly 
from those involving the sbottom vev. The remaining cubic terms, breaking supersymmetry 
softly, give constraints that are less tight than those obtained from UFB directions [10]: for 
this reason in this paper we focus on the latter constraints. 

In summary, the set of nonzero charged real vevs we work with is defined by: 

(h) = (in) = ^ , (h) = a(b R ) = i . (H/) = { H~ d ) = ^ . (T + > = (T~) = ^ . 

( 12 ) 


where we introduced a phase a = ±1 for later convenience, while the neutral ones are: 


<!/1> s/2 ’ ^ s/2 


(13) 


s/2 ’ ' ' s/2 

where we used a labeling for the neutral vevs different from that in Eqs. (9) to distinguish 
the vevs associated with the EW minimum from the unrealistic ones. Moreover, to simplify 
4 The stop vev turns out to be zero along UFB directions because of its Yukawa coupling, and therefore 


indeed only the sbottom labels need to be changed. 
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our following results without loss of generality, in the rest of the paper we assume v ho and 
all the Yukawa couplings to be positive. 

In the next subsection we determine which sets of vevs, among those defined in Eqs. (12,13), 
allow for a UFB direction in the potential. 


III.l. Relevant Vevs 


To determine the sets of nonzero vevs which allow for UFB directions, we look for those 
vevs combinations that can cancel all the D and quartic F terms. Assuming all the masses 
and couplings in Eq. (2) to be nonzero, requiring the superpotential derivative with respect 
to the triplet components in Eq. (4) to cancel, we obtain: 
cAUtessm 


= 0 , (t> C j = T° , T + , T => v H ± = 0 A (v H o = 0 V v H o = o) 


(14) 


Besides vh±, Eq. (14) requires either v H o or v H o to be zero. Indeed it can be shown that 
there is no vevs combination canceling all the quartic terms for nonzero v H o , and therefore 
we impose: 

v H o = 0. (15) 

Having defined v^o to be nonzero, we notice that m^ 0 , being large and negative at the EW 
scale for data points that feature a viable EW minimum, can generate a deep UFB direction. 

After imposing Eq. (15) and requiring the cancellation of the quartic F terms corre¬ 
sponding to the H° and fields, also the stop and charged triplet vevs turn out to be 
zero: 


cAUtessm 

dq 


= 0 , 


= H° u ,H d 


Vf = vt± = 0 . 


(16) 


Having set to zero the charged doublet and triplet Higgs vevs as well as the stop and the 
neutral down Higgs ones according to Eqs. (14,15,16), the only nonzero D and quartic F 
terms left are, respectively, 


Vo D S -^ (vf - vl 


s 2 \ f y 2 \ 

+ v h») 1 ^ D 4 (y 2 b v t + V2y b a\vlv H ov T o + — i4 g 4j . (17) 

Assuming a nonzero neutral up Higgs vev, we find therefore that it is possible to cancel all 
the quartic terms for the following sets of nonzero vevs: 

Vho ^ 0 A Vo 7 ^ 0 A ((vto = 0 A v~ h = 0) V (u^o = 0 A v~ b oc y/VH°) V (■ Vto ^ 0 A v~ h ^ 0)) , 

(18) 
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with the other vevs being all identically zero. The only other possible set of non-trivial vevs 
canceling all the quartic terms is 

Vto 7 ^ 0 V Vt± 7 ^ 0 , (19) 

with the other vevs, including u#o, being all identically zero. Evidently the potential in the 
UFB directions identified by Eq. (19) has the same form for either of the two vevs and one 
can work just with u^o. 

In the following subsection we work out the expressions for the potential along the four 
UFB directions expressed by Eqs. (18,19) and define the stability constraints associated with 
each of them. 

III.2. Stability Constraints 

We start with the simplest case, the UFB direction defined by Eq. (19): by setting to zero 
every vev but Vto in the potential (the corresponding result for Vt± is the same), Eq. (5), 
we get 

v 2 

fdjFB-i = \ m T + 2 Ht (B t + 2 fi T )\ . (20) 

We evaluate Eq. (20) at a renormalization scale A of the order of the heaviest mass in the 
physical particle spectrum [ 10 ], as to minimize the contribution of quantum corrections, 
which we neglect entirely. For the potential above the largest mass is roughly equal to v T o 
multiplied by the largest of its couplings, generally either or A. Rather than simply 
requiring the coefficient of the squared vev to be positive, a given point in parameter space 
is stable against tunneling 5 to UFB-1 if 

Vew ( v w) < Kjfb-i (A) , v w < A < A uv , v^o ~ 2 max \g 2 L , A 2 ] A 2 , ( 21 ) 

where all the couplings and dimensionful parameters are evaluated at A. For this purpose 
we calculated the full set of beta functions, including those of the dimensionful parameters 
which, to the best of our knowledge, were not given in the literature and that we present in 
Appendix A. In Section IV we elucidate how to implement in practice the stability constraint 
in Eq. (21) and those that follow in this Section. 


5 In this analysis we simplify the stability constraints by assuming the age of the universe to be infinite. 
The numerical factor in Eq. (21) agrees with [10], though for the numerical analysis it has little relevance. 



Next we take up the slightly more complicated case with only two nonzero vevs, the first 
one in Eq. (18): 


V\v H o^O,vz^O 2 


rnjvl + (m 2 H o + /^) v 2 H o + 


9y "k 9l ( 2 2 

-^6“ \ v * - v *2 


( 22 ) 


where m| is the soft mass squared of the slepton doublet L. By requiring the potential 
above to be flat along the v„ direction we find 


dv>V\ 


v„o 


= 0 


2 2 
v o = v m 


8 m 2 L 
9y + 9l 


(23) 


which, upon substitution in Eq. (22), gives the deepest UFB direction corresponding to our 
choice of nonzero vacua 


E 




UFB —2 — 


2 i 2,2 

"b + /i£) 


2 ml 
9y + 9l 


(24) 


Analogously to the stability constraint derived from the UFB-1 direction, a point in the 
TESSM parameter space may feature a stable EW minimum only if 


Uw M < VWb _2 (A) , Vw <A<A UV , v\ > 0 , t>J,j ~ 2 max [gl,X 2 , i/f] ‘ A 2 , 

(25) 

where all the couplings and dimensionful parameters are evaluated at A. 

The case with three nonzero vevs, the second one in Eq. (18), is a little more complicated, 
but the potential can be readily simplified by imposing its derivative with respect to Vp to 
be zero: 


d Vp V\ v 0 ^o,vr^o,vp^o 0 


2 2,2 
v o = + Vr - 


8 m 2 L 
9y + 9 2 l 


V\ 


V H o^=0,Vr^O,Vp^O 


= -T v t + -W ( m L + m Q + m l - V2y h avm9n) + ( m 2 L + m 2 m + y 2 D 


2m\ 
9y + 9l 


(26) 


By requiring the potential above to be flat along the v~ b direction and solving for v b we find 


«? = 


V2v h o y b \fi D \ - ml - m 2 Q - m? 


Vb 


(27) 


where, with our assumption that u#o and all the Yukawa couplings are positive, we simply 
took a equal to the sign of y D . In turn the vevs in Eqs. (26,27) identify the deepest UFB 
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direction corresponding to our choice of nonzero vacua 


Kjfb-3 — 


2 i 2 

m L + m Hl o 


2 m\ 


J H o 


+ 


\V 


D 


m 2 L + uiq + m~ 




V^yb 


ml + m 2 Q + mfj 

4 


.9y + f/i 


(28) 


where m 2 q is the soft mass squared of the squark doublet Q. The corresponding stability 
constraint that any point in the TESSM parameter space has to satisfy is 


Vew (««.)< V' UFB _3 (A) , v w < A < A uv , v?, > 0 , ~ 2 max [?/} , A 2 , y ' l ] 1 A 2 . 

(29) 

Notice that Vufb -3 is equal to the corresponding result in [10], plus the second from last 
term in Eq. (28), which is large and negative: this is because the authors in [10] determined 
vi by requiring the //[) quartic F term to cancel, rather than by solving a minimization 
condition. In case V}jo is small, v~ h turns out to be imaginary, and so we can require the 
potential in Eq. (26) to be flat along the v H a direction instead, which implies 


«f = 


V2v h o 


Vb 


m L + 7JI ho + l l n 


I^D 


(30) 


In this case the potential along the deepest UFB direction is 


V' 

v UFB-3 


mi + m 


H° 


m L + 


+ 1*d) 


2/1 


v m + 


D 


m L + m H° + l X D 


^D 


m 


L+ m Q + 


V^Vb 




9y + 9l 


(31) 


and the corresponding stability constraint reads the same as that in Eq. (29) but with Eu FB _ 3 
replacing Vufb- 3 - Notice that, contrarily to the vev in Eq. (27), the one in Eq. (30) is not 
generally smaller than v H o, and so one might be underestimating the heaviest mass, which 
is of the same order of the renormalization scale. For viable points, though, ml 0 is generally 
negative, in which case turns out to be of the order of v h[\ or smaller. For this reason we 
still determine Vffo according to the last one in Eqs. (29) and then vi by Eq. (27). 

Finally, we take up the last scenario in Eq. (18): for v ho ^ 0 A ^ 0 A Vto ^ 0 A v~ h ^ 0, 
the requirement for the potential to be flat along the v direction determines Vo as given by 
Eq. (26). The remaining minimization conditions produce rather involved solutions, which 
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turn out to be complex on a large and disconnected region of field space. For this reason we 
choose simply to cancel the quartic F terms, which is achieved by setting 


vl = 


\\v T a\v H o 

vVa 


(32) 


The potential along the plane identified by Eqs. (26,32) expressed in terms of the remain¬ 
ing two vevs is 


Fufb-4 — [m-T + 2 iir ( Bt + 2/ir)] 4— (jn? L + m^o + y 2 D j 

|Auj'| UjyO , 2 2 , 2\ 

+ _ ^ \rn L + m Q + m- b ) . 


2m\ 
9y + 9l 


2\/2y h v "^ ' bJ ' ^ 

By comparing the Erst line in the equation above with the Erst two UFB directions, 
Eqs. (20,24), and then realizing that the term in the second line is generally positive for 
phenomenologically viable data points, it is clear that any viable data point satisfying the 
first two stability constraints, Eqs. (21,25), is stable against tunneling to a vacuum along 
the plane defined by Eq. (33). 6 For this reason in this analysis we disregard the UFB-4 
stability constraint. 

In the next section we test the stability of TESSM at data points that satisfy the most 
relevant phenomenological constraints from experiment. 


IV. PHENOMENOLOGICALLY VIABLE PARAMETER SPACE 


The first relevant phenomenological constraint is given by the T parameter [23], which 
in TESSM receives a nonzero contribution already at tree level [L, 2] 


a P T = 


Srril 


w 


4v? 


m 


< 0.2 , 


2 2,2 

V = v u + v d , 


vi = v 2 + 4vl = 246 2 GeV 2 


(34) 


w 


where a e is the fine structure constant, the experimental constraint is at 95%CL [24], and 
the vevs appearing in the expression above are those defined in Eqs. (7). To satisfy the 
constraint above we take a small but non-zero fixed value for vt 


vt = 3\/2 GeV . 


(35) 


6 Notice that the last term in Eq. (33) is not allowed to turn large and negative upon renormalization for 
points featuring a viable EW vacuum, because otherwise the stability constraint in Eq. (29) would not be 
satisfied. 
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We then scan a large region of the TESSM parameter space, defined by 


1 < tp < 10 , 5 GeV < \fjbjj, ht\ < 2 TeV , 50 GeV < \M 1 ,M 2 \ < 1 TeV , 
| A t , A Tl B d , B t | < 2 TeV , 500 GeV < tuq, m f ~, < 2 TeV , 


(36) 


for data points producing the observed mass spectrum for SM fermions and gauge bosons 
and satisfying the direct search constraints defined below 



(37) 


where the neutral scalar masses are calculated by using the TESSM complete one loop 


effective potential [4], while the others are tree level masses. The above bounds on the mass 


of pseudoscalars and neutralinos are actually tighter than the experimental ones [24], as to 
avoid the phenomenological complications of invisible decays of the light Higgs, which are 


though relevant for dark matter [25]. For each of the scanned data points the couplings 
retain perturbativity at least up to 10 4 TeV, where we define a coupling constant to be 
perturbative if its renormalized value (obtained by running the given coupling by two loop 
beta functions [ ]) is smaller than 2n. ' Moreover the experimental constraints on a heavy 


neutral Higgs mass is imposed by rescaling the CMS bound on the mass of a SM heavy 
Higgs decaying to ZZ, m H o > 770 GeV [27]. 

After collecting 10957 viable data points satisfying the constraints outlined above, we 
test further their viability by imposing the stability constraints derived in the Section III. 

To perform the stability test we pick 100 energies equally spaced on a logarithmic scale 
between the EW scale, 246 GeV, and the chosen UV scale, A uv = 10 4 TeV, and run all 
the couplings and dimensional parameters to the 100 renormalization scales, with initial 
conditions defined by the given viable data point in the TESSM parameter space. For each 
of the 10957 viable data points we then evaluate the potential at the EW vacuum, Eq. (9), 
as well as at each of the 100 values along each of the UFB-1,2,3 directions, plus the one 
obtained from UFB-3, by replacing sbottom masses and couplings with the corresponding 
slepton quantities, which we call UFB-3'. The (real) vevs are automatically determined 
as functions of the given scale for each UFB direction. It turns out that only 24% of 

7 We choose 27 t because at larger values the couplings reach a fixed point which is an artifact of the series 
being truncated after the two loop contribution. 
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the 10957 viable data points satisfy the stability constraints defined in Eqs. (21,25,29). By 
applying each stability constraint individually one can assess the tightness of the constraints 
relative to each other: 58% of points are stable against tunnelling to vacua along the UFB-1 
direction, 95% against UFB-2, 67% against UFB-3, and just 41% of the scanned data points 
feature an EW minimum deeper than vacua along the UFB-3' direction, for scales up to 
10 4 TeV. Evidently the UFB-3' constraint, defined by Eq. (29) with slepton masses and 
couplings replacing the corresponding sbottom quantities, is by far the tightest among the 
3+1 constraints we imposed. This result is analogous to the one obtained for MSSM [10]. 
We also calculated the one loop potential evaluated at the EW minimum, and found it to 
be shallower than the tree level potential in Eq. (9) for 76% of data points, which makes the 
stability constraints defined by Eqs. (21,25,29) generally conservative. 

In Figs. (1,2) we plot all the 10957 scanned data points projected on two slices of the 
TESSM parameter space, m 2 Hd — mj lu plane (Fig. 1 left panel) and /i£> — Be plane (Fig. 
2 left panel), with grey points being unstable against at least one among the UFB-1,2,3,3' 
directions, while the colored ones satisfy all the stability constraints, Eqs. (21,25,29). The 
color code, shown in Fig. (1) (right panel), represents the value of the triplet coupling A 
at the given data point. Besides the requirements for an unstable trivial vacuum state and 



FIG. 1. Scanned TESSM data points satisfying all phenomenological constraints elucidated in the 
Section IV, projected on the m# — m 2 Hu plane (left panel). The grey points do not satisfy at 
least one of the stability constraints in Eqs. (21,25,29), while the viable data points are colored 
according to the value of the triplet coupling A as shown in the right panel. 

positive scalar squared masses, which make large regions of parameter space inaccessible, 
the stability constraints evidently rule out the region featuring a large negative soft squared 


13 





mass m 2 Hd as well as a large \hd\, with both parameters in absolute terms being smaller than 
1 TeV 2 and TeV, respectively, for stable data points. This in turn limits the mass of the 
lightest chargino and neutralino. 

In Fig. (2), right panel, we plot the TESSM cross section relative to its SM value for a 
light Higgs boson, produced at LHC and then decaying into a diphoton, as a function of the 
lightest chargino mass. Also shown are the average value (solid line) measured by ATLAS 
and CMS [28, 29], and the lower la bound (dashed line). The ratio of cross sections for the 



H D (TeV) m xt (GeV) 


FIG. 2. Scanned TESSM data points satisfying all phenomenological constraints elucidated in the 
Section IV, projected on the hd — Bo plane (right panel), and Higgs decay rate to diphoton as 
a function of the mass of the lightest chargino (right panel). The grey points do not satisfy at 
least one of the stability constraints in Eqs. (21,25,29), while the viable data points are colored 
according to the value of the triplet coupling A as shown in Fig. (1) (right panel). 


lightest Higgs boson, produced at LHC, decaying into ij particles is defined by 




0tot Hr ^ 





fl=h,qqh,... 


(38) 


where Br^ is the lightest Higgs boson branching ratio into the ij particles, a^ the production 
cross section of the given final state fl, and en is the corresponding efficiency. The details of 
the calculation of the branching ratio to diphoton and Higgs production total cross section 
in TESSM are given in [26], and we do not repeat them here. From Fig. (2) it is apparent 
that the stability constraints require the mass of the lightest chargino (and neutralino) to be 
lighter than about 700 GeV, for the scanned data points. As a consequence, within TESSM, 
a large deviation from the SM prediction is likely to be observed. 

Since a lighter mass spectrum is usually associated in supersymmetric theories with lower 
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fine tuning (FT), it is interesting to look also at the effect of the stability constraints on 
the heavier stop mass together with the fine tuning at each scanned data point. A simple 
estimate of FT in supersymmetry is given by the logarithmic derivative of the EW vev v w 
with respect to a given model parameter /x p [30, 31]: this represents the change of v w for a 
100% change in the given parameter, as defined below: 


= dlogvj 
^ ~ d log (A) 


, f4 (A) = nl {M z ) + log 


Pn 

167r 2 


A 

Mz 


^ = 16 AS’ (39) 


where in parenthesis is the renormalisation scale of (i p . The FT in m# , defined to be equal 
to f m H , is then given by [26] : 
log (A /M z ) 


FT = 


167T(9 ? ,2 m 2 rr 

u w Hu 


(6y 2 A 2 + 3A 2 A^ + 3\ 2 rn 2 Hd + 3A 2 m% + 3A 2 m 2 Hu — 2 g Y M 2 — 6g 2 L M,; 


+ 6 m 2 Q y 2 + Grrify 2 + 6 m 2 Hu y* + g 2 Y (3m? - m 2 Hd - 3m| + 3 m 2 Q - 6m? + m 2 Hu + 3m?)) , 


(40) 



where the derivative in the denominator acts on the expression of m 2 H , Eqs. (8). In Fig. (3) 


1.4 
m }2 (TeV) 


FIG. 3. FT as a function of the heavier stop mass (left panel) and tan/3 (right panel). The grey 
points do not satisfy at least one of the stability constraints in Eqs. (21,25,29), while the viable 
data points are colored according to the value of the triplet coupling A as shown in Fig. (1) (right 
panel). 

we plot the scanned data points as a function of the heavier stop mass (left panel) and 
tan/3 (right panel). The grey points do not satisfy at least one of the stability constraints 
in Eqs. (21,25,29), while the viable data points are colored according to the value of the 
triplet coupling A as shown in Fig. (1) (right panel). An interesting (and unwelcome) effect 
of the stability constraints is that the points with the least amount of FT are actually ruled 
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out, and as a result FT is on average 26% higher after imposing the stability constraints, 
somewhat offsetting the advantage for naturality of the tree level triplet contribution to 
the MSSM Higgs boson mass. With increasing heavier stop mass the amount of FT on 
average increases as well, as expected, and a larger portion of data points satisfy the stability 
constraints. Indeed also for small values of tan/3, or equivalently large values of the triplet 
coupling A, data points featuring a very low amount of FT are unstable against tunneling 
from the EW potential minimum to a vacuum along (at least one of) the UFB-1,2,3,3' 
directions. 


V. CONCLUSIONS 

In this paper we studied the Unbounded From Below (UFB) directions in the potential 
of a Y — 0 SU(2) triplet chiral superfield extension of MSSM, also called Triplet Extended 
Supersymmetric Standard Model (TESSM), and the associated stability constraints on the 
model’s parameter space. After introducing the model, we systematically looked for sets 
of nonzero vevs that can cancel the quartic terms belonging to the D and F sectors of the 
TESSM tree level potential, under some rather general and reasonable simplifying assump¬ 
tions, like chiral and charge symmetry of the vevs. We found four inequivalent sets of vevs 
that allow for UFB directions in the tree level potential, which we labeled UFB-1,2,3,4, 
respectively. One more UFB direction (UFB-3 7 ) is obtained by simply switching the sbot- 
tom couplings and masses in the potential, defined along UFB-3, with the corresponding 
stau quantities. Among those directions, UFB-2,3,3' turn out to be entirely equivalent to 
those already found in the MSSM potential [ LO] , in the sense that they do not contain any 
triplet contribution. For viable points in the TESSM parameter space, moreover, UFB-4 
is lifted if UFB-1,2,3 are lifted as well, which makes the UFB-4 stability constraint irrele¬ 
vant. The relevant stability constraints require the tree level potential evaluated at the EW 
minimum to be deeper than at any point along UFB-1,2,3,4,3', with couplings and masses 
renormalized at a suitable scale minimizing the (neglected) one loop contributions. To carry 
out renormalization, we furthermore derived the one loop beta functions for the TESSM 
dimensional parameters, which were not given in past literature. 

To assess the relevance of the stability constraints for TESSM, we first scanned the 
TESSM parameter space and collected 10957 data points which produce the observed SM 
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mass spectrum and satisfy direct search constraints on their superpartners and heavy Hig¬ 
gses, as well as EW precision parameter and perturbativity constraints. We then further 
tested the viability of these data points by checking how many of them satisfy the UFB- 
1,2,3,3' stability constraints: among the (otherwise viable) 10957 data points, only 24% 
turned out to be actually stable, with UFB-3' giving the tighter constraint, which rules out 
59% of the scanned data points. Two of the parameters that get constrained the most by 
stability are the soft up (Higgs mass, m# d , and the Higgs doublets supersymmetric mass, 
He >, both generally smaller than about 1 TeV for the data points featuring a stable EW 
minimum. As a consequence both the lightest chargino and neutralino turn out to be lighter 
than about 700 GeV for the same set of viable data points. While these are eventually 
too heavy for detection at LHC, one observable effect is that TESSM stable points feature 
on average a larger deviation from the SM predicted Higgs decay rate to diphoton than 
unstable points. We have shown furthermore that fine tuning within TESSM is on average 
26% higher after imposing the stability constraints, somewhat offsetting the advantage for 
naturality of the tree level triplet contribution to the MSSM Higgs boson mass. 

Finally, we conclude by highlighting the fact that the stability constraints described in 
this paper, by ruling out a large portion of the model’s parameter space and affecting the 
superpartners mass spectrum, should be taken into account in any phenomenological study 
of the triplet extension of MSSM. 


ACKNOWLEDGEMENTS 

This work was partially supported by the Academy of Finland, project 267842. MD 
acknowledges the hospitality of the Department of Theoretical Physics, Indian Association 
for the Cultivation of Science during the course of this work. SR acknowledges the hospitality 
of the University of Helsinki and the Helsinki Institute of Physics during the final stages of 
this work. 


17 



Appendix A: Beta Functions 


The beta functions at two loops in TESSM for the Yukawa couplings y t , y b , y T , A, as 
well as for the gauge couplings g 3 , g 2 = gL,gi = \/5/3 gy were already given in [ 6 ], using 
the same superpotential and soft terms as in Eqs. (2,6), so there is no need to write them 
again here. The dimensionful couplings are yx, y h , Mi, M 2 , M 3 , bx, bn, hr, h t , hb, h T , as 
well as the squared soft mass parameters m;L m 2 ~ , rn 2 - , m 2 -, m\ , m], rn 2 f , rrij . Their beta 
functions at one loop are defined by 


dz x 

dt 


a 

for = nr, y H , M u M 2 , M 3 , b Tl b H , h T , h u h b , h T ; t = log 


167T 2 


A 


EW 


(Al) 


and by 


dz x 

dt 


A 


(i) 


167T 2 


n 22222222 

for m fi , m d , m L , m~ e . 


(A 2 ) 


In the renormalization scheme using dimensional reduction (see [32] and references therein) 
with modified minimal subtraction (DR) we find 


Phi ~ ~ — 7g 2 hr + 3 fiTUb + ShrUb + dry 2 + -gfMiXx + 14<7 2 M 2 At + Lh b y b Xx + Qh t y t Xx 


2 . 6 2 


Aa? = -Hg'iht ~ 3g 2 h t - ^glht + h t y£ + —g 2 M\y t + 6g 2 M 2 y t + ^g^M 3 y t + 2h b y b y t + 18 h t y 2 


+ 2h r y r Xx T 24/i-j’Aj' , 

13 2 1 0„2 7 16 2 


(A3) 


26 


32 


ht is* 1 "" 3 

6h T y t X T + 3Mr 7 
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glK - 3g 2 2 h b - + 6g 2 M 2 y b + M 3 j/ b + 18/^ + 2hty b y t + 


(A4) 
2 


+ 2 h T y b y T + + 6h T y b X T + 3h b X 2 T , 


(AS) 


m 9 18 

= -pi\K - 3g\h r + 3/i r ?/ 2 + —g 2 l M l y T + Gg 2 2 M 2 y T + 6 h b y b y T + 12 /r T ?/ 2 + 6h T y T X T + 3 h T Q$6) 


A$ = dr (Sgl + 4 A^) , = dH (-^1 - 3 ^2 + 3 2 /& + % 2 + vl + 6 A|^ , 


A$ = , /3g = 6^ 2 2 M 2 , ^ = -6glM :i , 

A&y = —8g 2 bx + 4:bxX^ + 16 g 2 M 2 yx + 8hxXxyx , 

?(i) 3 


(AT) 

(AS) 

(A9) 


6 


Pb’ = ~^glb H - 3g]b H + 6h b y b dH + 3b H y 2 t + 6b H Xl + b H y 2 + 3b H y 2 b + -g^Mxyn + 6g 2 2 M 2 y H 


Sht/myt + 12hx^HXx + 2h T yHy T ■ 


(A 10 ) 


By defining the quantity 


S = m\ — m 2 f, + 3 m 2 A — 6rri 2 + 3m] — 3m] + 3rti 2 , 

M u Md Li a a L e 


(All) 
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we can write the square mass parameters’ beta functions as 


$> = 


4Aji?n'h — 16y 2 Af 2 4:hx^T^T ~b 4/4 4- 4A j^ui 2 £ 
3 « 


l t a t 

Pfil = \v\ S 6 7 8 9 10 11 12 + 6A T™| d - 4 ^.9 i M i 2 - 6y 2 M 2 2 


2m|,Ay , 


+ 4m|A^ , (A12) 

' 2 Af 2 2 + 6 h 2 + 6/4 + §yl m \ u + 6A^m| u + Gmjy 2 + 6 m?y 2 


(A13) 

4 m| d + 2 y 2 m| d + 2 m 2 y 2 - jj 

2/r 2 + 6/4 + 6A| 


/?2! = % 2, "4 d + %6 m l + 6/ ^ + 6 ^ m 0 + 6 A t"4 


^ 1} = 2yf™ 2 




4 1 ’ 


(!) 

d 

(1) 

L 

(1) 


4m^ + 2m|y 2 + 2m 2 f \ 2 T 
ib rr > 2 H d + 2y 2 m 2 ~ + 2 h 2 b + 2y 2 m| - 2 
2-m|y 2 + 2m 2 y 2 , 

-^i M i 2 _ y^ 3 M 3 - ^5'i‘S' + 4/r 2 + 4y 2 m 2 fu + 4mjy 2 + 4m 2 y 2 , 
%&™| d + 4y 2 m 2 + Ah 2 b + 4y 2 m| - y y,M 2 - y y|M 3 2 + , 

2 y 2 r m \ d + 2m 2 y 2 - y 2 M 2 - 6g 2 2 M 2 - y 2 Y + 2 h 2 r + 2m|y 2 , 
4y 2 m 2 > + Amjyl - ^o 2 M 2 + ^g?S' + 4h 2 + 4m 2 f y 2 . 


ly 2 - - 6 g 2 M 2 - | 

yyf M 3 2 + 5 5 'i‘^ + 2h 2 


5^ 

(A14) 

b 2 + 2y t 2 m^ 
(A15) 

(A16) 
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